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ABSTRACT 

Does the "blazar sequence" exist, or is it a result of a selection effect, due to the difficulty in 
measuring the redshifts of blazars with both high synchrotron peak frequencies (> 10 15 Hz) and 
luminosities (> 10 46 erg s^ 1 )? We explore this question with a sample of blazars from the Second 
Catalog of Active Galactic Nuclei (AGN) from the Fermi Large Area Telescope (LAT). The Compton 
dominance, the ratio of the peak of the Compton to the synchrotron peak luminosities, is essentially a 
rcdshift-independent quantity, and thus crucial to answering this question. We find that a correlation 
exists between Compton dominance and the peak frequency of the synchrotron component for all 
blazars in the sample, including ones with unknown redshift. We then construct a simple model 
to explain the blazar properties in our sample, where the difference between sources is due to only 
the magnetic field of the blazar jet emitting region, the external radiation field energy density, and 
the jet angle to the line of sight, with the magnetic field strength and external energy density being 
correlated. This model can reproduce the trends of the blazars in the sample, and predicts blazars 
may be discovered in the future with high synchrotron peak frequencies and luminosities. At the same 
time the simple model reproduces the lack of high-synchrotron peaked blazars with high Compton 
dominances (> 1). 

Subject headings: galaxies: active — BL Lacertae objects: general — quasars: general — gamma rays: 
galaxies — radiation mechanisms: nonthermal 



1. INTRODUCTION 

AGN with relativistic jets pointed along our line 
of sight are known collectively as blazars (e.g., 
lUrrv fc PadovanH fl995). Blazars have two main sub- 
classes, those with strong broad emission lines, Flat 
Spectrum Radio Quasars (FSRQs), and those with weak 
or ab sent lines, known a s BL Lacertae ob jects (BL 
Lacs; iMarcha et all 119961 lLandt et all 120041). FSRQs 
are t hought to be Fanaroff-Riley (FR; iFanaroff fc Rileyl 
119741 ) type II radio galaxies aligned along our line of 
sight, while BL Lac objects are thought to be aligned 
FR type I radio galaxies. In general, the spectral energy 
distributions (SEDs) of blazars have two basic compo- 
nents: a low frequency component, peaking in the optical 
through X-rays, from synchrotron emission; and a high 
frequency component, peaking in the 7 rays, probably 
originating from Compton scattering of some seed pho- 
ton source, either internal (synchrotron self-Compton or 
SSC) or external to the jet (external Compton or EC). 

Aside from their clas sifications as FSRQs or BL Lacs 
from optical spectra, lAbdo et"aT1 (|2010cD subdivided 
them based on their synchrotron peak. They are consid- 
ered high synchrotron-peaked (HSP) blazars if their syn- 
chrotron peak v pk > 10 15 Hz; intermediate synchrotron- 
peaked (ISP) blazars if 10 14 Hz < v s p v k < 10 15 Hz; and 
low synchrotron-peaked (LSP) blazars if v s pk < 10 Hz. 
A lmost all FSRQs are LSP blazars. 

iFossati et alj (|1998t ) combined several blazar surveys 
and noted an anti-correlation between the luminosity at 
the synchrotron peak, L s p v kl and the frequency of this 
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peak, v pk . They also noticed anti-correlations between 
the 5 GHz luminosity (L5 GHz) and v pk ; the 7-ray lu- 
minosity and v pk ; and the 7-ray dominance (the ratio of 
the EGRET 7-ray luminosity to the synchrotron peak 
luminosity) and v s pk . These correlations were claimed 
as evidence for a "blazar sequence" , a systematic trend 
from luminous, low-peaked, 7-ray dominant sources with 
strong broad emission lines to less luminous, high-peaked 
sources with weak or none xistent broad emission lines 
and 7-ray dominance ~ 1. Ghiscllini et al. (1998) pro- 
vided a physical explanation for these correlations. If 
the seed photon source for external Compton scattering 
is the broad- line region (BLR), and the BLR strength 
is correlated with the power injected into electrons in 
the jet, one would expect that more luminous jets have 
stronger broad emission lines and greater Compton cool- 
ing, and thus a lower v s pk . As the power injected in 
electrons is reduced, the broad line luminosity decreases, 
there are fewer seed photons for Compton scattering, and 
consequently the peak synchrotron frequency moves to 
higher frequencies. This trend is also reflected in the 
lower luminosity of the Compton-scattered component 
relative to the synchrotron component as v pk moves to 
higher frequencies. If this physical explanation is cor- 
rect, it provides a powerful tool for understanding blazars 
and their evolution, not unlike the "main sequence" for 
stars on the Hertzsprung- Russel diagram. The blazar se- 
quence also has implications for "feedback", a relation- 
ship where the AGN jet heats the hot, X-ray emitting 
intracluster medium (ICM), while the ICM gas Bondi ac- 
cretes icjito_tJ-ieJ^ck_^ jet 
(e.g.. IBirzan et al J 120041 ). IHardcastle et al.1 (|2007fl have 
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suggested that radio galaxies without high excitation 
narrow lines (mostly FR Is, and presumably BL Lacs) 
are fed by "hot mode" accretion, i.e., accretion from the 
hot ICM; while radio galaxies with high excitation nar- 
row lines (mostly FR lis, and presumably FSRQs) are 
fed by "cold mode" accretion, i.e., accretion of cold gas 
unrelated to the ICM. 

The L 8 pk -v s pk anti-correlation has been questioned. Us- 
ing blazars from the Deep X-ray Radio Blazar Survey 
and the ROSAT All- Sky Survey/Green Bank Survey, 
iPadovani et al.l ([20031 ) did not find any anti-correlation 
between v s pk and radio, BLR, or jet power. This 
work, however, has been criticized for its relatively poo r 
SEP cha r acter ization (Ghisellini fc Tavecchiol [2008b). 
iPadovariH ([2007ft has identified several major predictions 
of the blazar sequence: the anti-correlation will continue 
to be found in more complete samples; since low lumi- 
nosity objects are almost always more plentiful than high 
luminosity objects, high-peaked blazars should be more 
plentiful than low-peaked blazars; and the lack of plenti- 
ful outliers, i.e., objects that are low-peaked and faint, or 
high-peaked and bright. If any of these predictions are 
c ontradicted, it would invalidate the sequence. 

iGhisellini fc Tavecchiol ((2008b) pointed out that since 
blazars are anisotropic emitters, the lower left of the 
diagram (i.e., sources that have low peaked, faint syn- 
chrotron components) should be filled in by sources 
which are viewed increasingly off-axis. But the lack 
of sources in the upper right region of the F s pk -v s pk 

plot (i.e., sources that have high-peaked, bright syn- 
chrotron co mponents) could be the result of a selec- 
tion effect (IGiommi et al.1 120021: IPadovani etall 120021 : 
iGiommi et al.l 20051 I20i2al) . Since a large fraction of 
BL Lac objects have entirely featureless optical spec- 
tra, their redshifts, z, and hence luminosities, are im- 
possible to determine. These could be extremely lu- 
minous, distant BL Lacs that would fill in the up- 
per right region. iNieppola etUI (|2006l ) found no cor- 
relation between the frequency and luminosity of the 
synchrotron peaks for objects in the M etsahovi Ra- 
dio O bservatory BL Lacertae sample. iChen fc Bail 
([20T1 did find an anti-correlation, using source s found 
in the LAT bright AGN sample (jAbdo et al.l l2009aL 
l2010cf K In recent works, an "L" -shape in the L^-v^ 
plot seems to have emerged, as lower luminosity low- 
peaked sourc es have been detec t ed with more sensitiv e 
instruments (iMever et al.l 120111: IGiommi et"all l2012bl ). 
INieppola et al.l ( |2006f ) found more of a "V" shape, al- 
though their plot did not include FSRQs; if high lumi- 
nosity and low-peaked FSRQs were added, it might ap- 
pear as more of an "L". None of the studies that find 
few sources with high L s pk and high v pk have yet to ad- 
dress the sources without redshifts, however. Rau et al.l 
<|2012f ) provide reliable photometric redshifts for 8 BL 
Lacs at z > 1.3, and four of these do inde ed seem to have 
v2 > 10 15 Hz and L s p y k > 10 46 erg s" 1 (|Padovani et al.1 

l20ll . 

With the advent of the Fermi Gamma-Ray Space Tele- 
scope era, it is now possible to also characterize the 
Compton peak frequency v pk and luminosity L pk for a 
larger number of objects than previously possible with 
EGRET. As we show, the Compton dominance Ac = 



L pk /L s pk is an important parameter for characterizing 
this sequence, since it is a redshift-independent quan- 
tity. We explore the sequence using a sample based 
on the second catalog of AGN from the LAT (2LAC; 
Ackcrman n et al.l 1201 ID in Section [2j We then show in 
Section [3] that the sequence, including sources with high 
L s pk and high v s pk can be reproduced with a simple model 
involving power-law injection of electrons and radiative 
cooling. Finally, we conclude with a discussion of these 
results (Section 2J . 

2. THE 2LAC BLAZAR SEQUENCE 
2.1. Sample Definition and SED Characterization 

The 2LAC (|Ackermann et al.l 1201 ID presents the 
largest 7-ray catalog of blazars yet. It allows for the 
characterization of the high energy component for a 
greater number of blazars than previously possible. Here 
we look at the blazar sequence among the 2LAC clean 
sample, which includes 886 total sources, with 395 BL 
Lacs, 310 FSRQs , and 157 sources of unknown type. 
lAckermann et all ([20111 ) collected the fluxes of many of 
these sources at 5 GHz, 5000 A, and 1 keV, and used em- 
pirical relations for finding the peak frequency of the syn- 
chrotron component from the slope between the 5 GHz 
and 5000 A flux (a rn ) , and between the 5 00 A and 1 
keV flu x (a ox ) found bv lAbdo et all (j2010d ). lAbdo et al.1 
(|2010d) in turn had found these empirical relations by 
fitting the broadband SEDs of 48 of the bl azars in the 
3-mon th LAT bright AGN sample (LB AS; lAbdo et al.1 
12009a!) with third degree polynomials to determine the 
peak synchrotron frequency, v p k . With the v s pk est imated 
from the empirical relations, lAckermann et all ([201 ID 
could classify the objects in the 2LAC sample as LSP, 
ISP, or HSP. Here we use their results for v?*, which 
are corrected for redshift, and so are in the sources' rest 
frames. The 2LAC also includes the 5 GHz flux den- 
sity for many sources. We combine this with a spectro- 
scopic redshift mea surement, i f availa ble, or photometric 
redshifts found by iRau et all ([2012), to get the lumi- 
nosity distance, <iJ3 and calculate the 5 GHz luminosity 
L5 GHz — ^d 2 L (i>F v ) 5 GHz- We have corrected the ra- 
dio luminosities to how they would appear in the rest 
frames of the sources (i.e., fc-corrected them), assuming 
a r = 0, where the radio flux density is F v oc v~ ar . A plot 
of L 5 GHz versus v s pk for the objects in the 2LAC clean 
sample for which the catalog has a liste d v s pk l 5 GHz flux 
density, and z (or a z from Ra u et aT1l2012D is given in 
Figure [TJ This includes 352 sources, of which 145 are 
BL Lacs, 195 are FSRQs, and 12 are AGN of unknown 
optical spectral type (AGUs; i.e., unknown whether they 
a re FSRQs or BL La cs). 

lAbdo etall (|2010cl ) also provided an empirical formula 
for determining the flux at the synchrotron peak {F pk ) 

from the the 5 GHz flux density and v s pk . We use this 
empirical relation to determine the peak synchrotron flux 
for the 352 objects in our sample, and again combine 
it with oIl to calculate the peak synchrotron luminosity, 
L s pk . A resulting plot of L s pk versus v s pk is shown in Figure 

m 

1 To calculate d^, we use a flat ACDM cosmology with Ho = 
71 km s -1 Mpc -1 , Q m = 0.27, and Q A = 0.73. 
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Figure 1. Radio luminosity at 5 GHz versus peak synchrotron 
frequency for the 2LAC clean sample. Symbols are the same as in 
Figure [2] 
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Figure 3. The LAT spectral index (r 7 ) and the peak frequency of 
the Compton component, u^ k de termined from third-order polyno- 
mial fits from Abdo ct al. (2010c), plotted as circles. The empirical 
fit they determined is plotted as the line. 
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Figure 2. Peak synchrotron luminosity versus peak synchrotron 
frequency for objects in the 2LAC clean sample. Filled circles rep- 
resent FSRQs, empty circles represent BL Lacs, and filled squares 
represent objects which do not have an unambiguous classifica- 
tion. Dashed lines indicate the boundary between HSPs and ISPs 
(u sy k = 10 15 Hz) and between ISPs and LSPs (v'Y = 10 14 Hz). 



Additionally, lAbdo et all (|2010cf) fit the high energy- 
components of their 48 LBAS blazars with a third de- 
gree polynomial to determine the peak of the 7-ray com- 
ponent (presumably from Compton scattering). They 
found an empirical relation between v pk and the LAT 



7-ray spectral index, T 7 . Figure [3] shows T 



from their polynomial fits to the objects in their sample, 
and the empirical relation they found. The empirical fit 
seems reasonable in the range 1.6 < T 7 < 2.6, but ex- 
tending the relation outside this range is questionable. 
Approximately 10% of the 352 sources are also found in 
the 58-month Swift Burst Alert Telescope (BAT) cata- 
log. For these sources, we extrapolated their BAT and 
LAT power-laws and found where the extrapolations in- 
tersected. If they intersected within the range 195 keV to 
100 MeV (that is, in between the BAT and LAT band- 



passes) we used this location as 



''pk " 



Using the BAT 



spectrum in this way allows an improved estimation of 

2 http : //heasarc . nasa . gov/docs/swif t/results/bs58mon/ 



Vpk over t ne empirical relation from I Abdo et al.l (|2010cl) . 
particularly for those very soft sources, for which this em- 
pirical relation is untested. For other sources, the peak 
was determined fro m the LAT sp e ctral in dex and the em- 
pirical relation from I Abdo et al.l ()2010cf) . In principle, a 
similar technique could be used by combining the LAT 
spectra and very-high energy (VHE) spectra taken from 
atmospheric Cherenkov telescopes, in order to find v^ k 
for hard LAT sources. However, the VHE spectra tend 
to be highly variable, and are not integrated over a long 
period of time, and often are not simultaneous with the 
LAT spectra. Since the BAT and LAT spectra overlap 
in time and are integrated over long times, it seems more 
reasonable that they would represent a similar state. 

Once the location of v pk is known, either by using the 
empirical relation or from the BAT-LAT intersection, the 
flux at the peak, F^ k , can be estimated by extrapolating 
the LAT spectrum, and the corresponding luminosity can 
be calculated with iS. = ^d 2 L F^ k . To test the accuracy 



pk 



of this approximatio n, we use this meth od to estimate 
Lpk f° r the sample in lAbdo et al.1 (|2010d ). and compare 
it with the value found by their fits to 48 sources. The 
results can be found in Figure [H the agreement seems 
reasonable. Therefore, we used this technique to estimate 
Lp k for the 352 objects in our sample. With this, we can 

calculate the Compton dominance, Ac = ±j pk 



tC ifsy a 



plot of Ac versus i/fj[ is shown in Figure [5] Also note 
that Ac is independent of redshift (Ac — L^/L 8 ^ £=; 
Fpk/Fpk)- Thus we can plot in Figure [5] an additional 
170 sources from the 2LAC clean sample that have well- 
determined synchrotron bumps but do not have known 
rcdshifts. For these sources, the plotted vji is a lower 



pk 

since the redshifts are not known. However, u pk 



limit 

will be larger by only a factor (1 
few. 



z), i.e., a factor of a 



2.2. Luminosity and the Sequence 
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log (L pk / erg s ) (from fit) 

Figure 4. The peak Compton luminosity, L^ k , estimated from 

empirical relations as discussed in the text, plotted versus L^ k 

determined from the fits by Abdo ct al. (2010c). The line shows 
where the estimated and fit values would be equal. 
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Figure 5. Compton dominance (i.e., L^ k /L^ v k ) versus peak syn- 
chrotron frequency. Symbols are the same as in Figures [T] and [2] 
additionally, rightward-pointing triangles represent BL Lacs with 
unknown rcdshifts, for which v sy , is a lower limit. 



In their diagrams of L sy k v ersus v^f., INieppola et al.l 

(|20€6l) and lMever etHI (l20ll note a "V" or "L" shape. 
This shape is not at all obvious in our version of this plot 
in Figure [5J To determine the strength of the correlation 
between L s p v k and we have computed the Spearman 
(p) and Kendall (r) rank correlation coefficients and the 
probability of no correlation (PNC) calculated from each 
coefficient. The results can be found in Table [TJ The 
results from p and r are similar in all cases. The PNC is 
very small for the BL Lacs and FSRQs separately, and 
even lower for all the sources combined, where the prob- 
ability is essentially zero that ther e is not a correlation . 
The small PNC is consistent with IFossati et al.l (|1998l) . 
who find essentially a zero chance of no correlation using 
the r test. Note however that sources with unknown z 
are not included. Since these sources could have high 
Ljf., their exclusion could explain the anti-correlations 
for the whole sample and for the BL Lacs in general 
(jGiommi et al.ll2012al ). since they could fill in the upper 
right part of this diagram, as mentioned in Section [TJ In- 



deed, several sources with redshifts estimated from pho- 
tome t ry have been fou nd in this region (IPadovan i et al.l 
[20lllRau et~aII [2012K and are included in our sample. 
The objects with unknown z are almost certainly BL 
Lacs (i.e., they probably have small broad-line equivalent 
widths, making their rcdshifts difficult to determine), so 
it would not affect the hint (92.7% and 93.4% for the p 
and t coefficients, respectively) of anti-c orrelation found 
for FS RQs. It should also be noted that Nicppol a et all 
( 2006) did not find a correlation between L s p v k and u pk for 

BL Lacs using the p coefficient, possibly because they 
used log-parabola f unctio ns to characterize the SEDs, 
while IFossati et all (|1998f ) and lAbdo et al.l (|2010d) use 
third-degree polynomials. Since our determinations of 
v f. and L s pk are based on empirical relations related to 
these third-degree polynomial fits, our re sults could also 
be exp ected to be different from those of INieppola et all 
(I2006T). 

IFossati et al.l (|1998f) also found a significant anti- 
correlation between the 5 GHz luminosity, L§ GHz, and 
Vpfc, using the r test. Like them, we find the anti- 
correlation with L5 ghz versus v% (Figured]) to be much 
more significant than the one between L s p y k and fp k , as 
can be seen in Table [TJ This is true for both FSRQs and 
BL Lacs alone as well as together. If the anti-correlation 
is ex plained by the increas ing cooling at higher luminosi- 
ties (jGhiscllini et al.lll998f ). one might expect the corre- 
lation with L S p k to be more significant, since the emission 
at 5 GHz is thought to be from a different region of the 
jet than the emission at the peak, and the region emit- 
ting at the peak should be synchrotron self-absorbed at 
5 GHz. However, since v^j, was determined in part based 
on the 5 GHz flux, the strong correlation with L5 GHz 
is probably a re sult of this dependence. As pointed out 
bv iLister et all poTTI ). if Vp k increases, but its spectral 
shape and vF v peak flux do not change, the radio flux 
(or luminosity) will naturally decrease. This can easily 
explain this anti-correlation. However, what is unclear 
is whether a peak frequency derived from the radio flux 
should be interpreted as a cooling break, since these two 
should be from different regions and possibly indepen- 
dent. Determination of v s pk independent of low radio fre- 
quency emission should be a good, although technically 
challenging, way to test this. 



2.3. Compton Dominance and the Sequence 



Combining 
Fossati ct al 



their 



nance versus v„ 



results with 
made a plot 



EGRET 
of 7-ray 



data, 
domi- 



y pk . Using LAT data from the 2LAC, as 
described in Section |2~TI we make a similar plot (Figure 
0, although we use L pk instead of simply the 7-ray 
luminosity. A distinct "L" shape is seen in this figure. 
We have also computed the correlation coefficients p 
and t for Ac versus v pk , and the results can be found 
in Table [TJ First we do this only for sources with 
known z. There is no evidence for a correlation for the 
FSRQs alone, although the probability that there is no 
correlation for the BL Lacs alone is very small. For 
the combined sample of all sources with known z, there 
is essentially zero chance that there is n o correlation, 
similar to the L sy k versus v sy k correlation. IFossati et al I 
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Table 1 



Statistics of correlations involving V k - 



Sample 


P 


PNC(p) 


r 


PNC(r) 


T sy sy 

L v \ versus i/ p j» 


BL Lacs 




-0.19 


0.020 


-0.12 


0.035 


FSRQs 




-0.12 


0.073 


-0.088 


0.066 


All sources 


with known z 


-0.54 


4.2 x IO" 28 


-0.35 


< IO" 50 


L 5 GHz versus i/° y k 


BL Lacs 




-0.64 


5.4 x 10" is 


-0.45 


< io- bU 


FSRQs 




-0.36 


1.82 x 10" 7 


-0.25 


1.79 x 10~ 7 


All sources 


with known z 


-0.79 


< 10~ 50 


-0.58 


< 10" 50 


Ac versus u"" k 


BL Lacs 




-0.30 


2.7 x 10" 4 


-0.21 


1.3 x 10" 4 


FSRQs 




8.9 x IO" 3 


0.90 6 


3 x 10~ 3 


0.89 


All sources 


with known z 


-0.66 


9.8 x IO" 45 


-0.45 


< IO" 50 


All sources 


l a 


-0.66 


< io- 50 


-0.46 


< 10- 50 


All sources 


2 b 


-0.66 


< io- 50 


-0.46 


< 10~ 50 


All sources 


3 C 


-0.64 


< io- 50 


-0.44 


< 10~ 50 


a All sources 


including those with unknown 


2, assuming they 


are at 2 — 


0.00. 


All sources 


including those with unknown 


2, assuming they 


are at 2 — 


0.35. 


c All sources 


including those with unknown 


z, assuming they 


are at 2 = 


4.00. 



(1998) also find a low probability that there is no 
correlation using the r test, although our result is 
stronger. We also computed the coefficients for all 
sources, including the ones with unknown z, computing 
their v s ^ k assuming z = ("all sources 1" in Table 

[J). We find essentially no chance that the addition of 
sources with unknown z could ruin the correlation when 
these sources are included. The term v^j. will vary by a 
factor of a few due to redshift, so we also calculated the 
coefficients assuming all these sources with unknown z 
are at z = 0.35, the average of the BL Lacs with known 
z ("all sources 2"); and assuming these sources are at 
z = 4 ("all sources 3"), which is significantly higher 
than the maximum redshift of the entire sample (which 
is z — 3.1). In each case, there is essentially a 100% 
chance that an anti-correlation exists. Although the 
correlation between L a ^ k and i/f could exist only because 
of a redshift selection effect, a similar redshift selection 
effect cannot explain the correlation between Ac and 
v pk ■ The relationship between Ac and v s ^ k thus seems to 
have a physical origin. 

Although a similar quantity, 7-ray dominance, was 
discussed as part o f the original "blazar sequence" by 
iFossati et al.l (fl998). 7-ray dominance or Compton dom- 
inance has been mostly overlooked ever since, in favor 



of testing the correlation between L S H 



and 



sy 



Padovani et al.l 120031: 
20071 IMever et al.11 2011 



Nicppola et al 



IGiommi et al 



n k (e-g-, 
[200l IPadovanil 



I2012af ). although 



see lGiommi et al.l (|2012bl ) for a brief discussion. This ne 
gleet may be in part due to a lack of quality 7-ray data, 
a deficiency that has been corrected in the Fermi era. 

Errors and Outliers 



2.4. 



The errors on determining 



sy 
pk 



and have essen- 
tially been ignored in the past. Although for radio 
and optical observations measurement errors are small, 
clearly the peak will depend on how well-sampled the 
synchrotron bump is. Furthermore, X-ray measurements 
can have large errors (due in part to the necessity of as- 
suming a spectral form to convolve with the response 
matrix of an X-ray instrument), and here ignoring their 



errors could lead to significant errors on and L e Jj,. 
It is also possible that the peak location can depend 
on the fitting function used. It does not seem to be 
stand ardized; sometimes a log-parabola function is used 
(e.g., [Nicppol a et al.l |2006[) . sometimes a t hird degree 
polynomial is used (e.g.. I Fossati et al.lll998t lAbdo et"al~l 
I2010d : IMever et al.l 120111) and sometimes a physically- 
motivated synchrotr on/ Compton model is used (e.g., 
iPadovani et al.l l2003f) . although the third-order polyno- 
mial fit is the most common. 

When locating the Compton peak, the measurement 
errors can have an even greater effect, since the 7-ray 
error bars tend to be larger than at lower frequencies. 
This can lead to s ignificant errors in determining Ac- 
lAbdo et al.l (|2010d) have estimated the error between the 
approximate expressions and the polynomial fits to be 
about an order of magnitude. This estimate still neglects 
the measurement errors, which can have a large effect on 
results. 

The two sources with unknown redshifts in the up- 
per right quadrant of Figure [5] are 2FGL J0059.2- 
0151 (1RXS 005916.3-015030) and 2FGL J0912.5+2758 
(1RXS J091211.9+27595) with LAT spectral indices of 



L-v = 1.15 ± 0.36 and 



r 7 = 



1.20 ± 0.37, respectively. 



These are the two hardest sources in the 2LAC, and the 
sources with the largest error bars on their spectral in- 
dices; only one source in the 2LAC clean sample is fainter 
than these sources (2FGL J1023.6+2959). They are 
clearly outliers. Propagating the error on their spectral 
indices, one finds that they have Compton dominances 
of log 10 (A c ) = 1.44 ± 1.82 and log 10 (A c ) = 1.48 ± 1.98, 
respectively; they have Ac consistent with unity within 
their error bars, and so are consistent with the "L" shape 
seen in Figure [5] 

There are several FSRQ outliers in Figure [5] as well; 
perhaps most interesting are those with Ac > 10 2 . 
These three sources are: 2FGL J1017. 0+3531 associ- 
ated with B2 1015+35B; 2FGL J1154.4+6019 associ- 
ated with CRATES Jl 154+6022; and the most extreme, 
2FGL J1522.0+4348 associated with B3 1520+437, with 
Ac = 1500, which is not shown in Figure [5] The 
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sources 2FGL J1017.0+3531 and 2FGL J1522.0+4348 
are estimated to have Compton components peaking 
in the BAT bandpass, both with peak fluxes > 6 x 
10~ 10 erg s _1 cm -2 , and both would then be visible 
with BAT. The source 2FGL J1154.4+6019 is estimated 
to have a high-energy peak observed at 577 keV. Ex- 
trapolating from this peak back into the BAT band- 
pass with a spectral index T — 1.5 gives a flux of 
9.6 x 10~ n erg s _1 cm~ 2 at 100 keV. There are several 
objects in the BAT catalog with fluxes less than this, so 
the BAT would probably have detected this source as 
well, if this estimate is correct. Thus, these three sources 
probably have a low-energy cutoff in their Compton com- 
ponent before the BAT waveband, and their Ac shown 
here are probably not accurate. 

Variability is another factor that can affect correlations 
determined here. Blazars are highly variable at all wave- 
lengths, and the p eaks of blazars can vary by several or- 
ders of magnitude dFossati et alj[2000t iCostamante et al] 
120011) . iGiommi et al.l (|2012b|) explore the issue of charac- 
terizing peaks at different epochs from microwave, X-ray, 
and 7-ray data. They conclude that microwave variations 
are rather minimal and make little difference in charac- 
terizing SEDs, but variations in X-rays (0.1 - 2.4 keV) 
and 7-rays (in the LAT energy band) can vary by up to 
a factor of 10 or greater. 

For all of these reasons, the results in Figures [U [2l 
and [5] are probably not accurate for individual sources, 
although no one has rigorously estimated their errors. 
However, we are more interested in the overall trend than 
individual sources, so the trend could still be present 
even for large errors. To test if it will still be present, 
we performed a simple Monte Carlo simulation. For 
each source, we randomly drew the LAT spectral in- 
dex (r 7 ) and luminosity (£ 7 ) from the measured val- 
ues and measurement errors, assuming the errors are 
described by normal distributions. Then we used this 
randomly-drawn value to calculate Ac- There are no re- 
ported errors on v^_, so we assumed that log(f^) has 
normally-distributed errors of one decade. We then used 
the Spearman and Kendall tests on the randomly-drawn 
distribution to determine if the correlation persists. We 
performed this simulation 10 5 times, and found that a 
PNC > 5.7 x 10~ 7 (the probability corresponding to 5<r, 
assuming normally distributed errors) was found in only 
around 5% of the simulations, for both p and r tests. 
Thus, we are confident that the trend persists despite 
large (random) errors. There do not seem to be any sys- 
tematic errors in using the empirical re lations instead 
of the polynomial fits (jAbdo et al.ll20l6"cl Figures [3] and 
0]), except possibly a small systematic underestimation 
of Lp k as shown in Figure |U 

3. THEORETICAL BLAZAR SEQUENCE 
3.1. Simple Unified Blazar Emission Model 

We describe a simple model for blazar jet emission, and 
show that it can reproduce blazar properties described in 
Section [21 This mod e l is sim ilar t o the ones presented by 
iBottcheTfe Dermerl (12002|) and iGhisellini fc Tavecchiol 
(2008b[): IGhisellini et all (j2009D . We assume the rela- 
tivistic jet is dominated by emission from a single zone 
which is spherical with radius R b in its comoving frame, 
and moving with highly relativistic speed /3c giving it 



a Lorentz factor T = (1 — /3 2 ) -1 / 2 . The jet makes an 
angle to the line of sight 9 so that the Doppler factor 
Su = [r(l — f3 cosO)}^ 1 . Electrons are injected with a 
power-law distribution given by 



Qe(7) = Qo7~* #(7571,72) 



(1) 



where the Heaviside function H (x; a, b) = 1 for a < x < b 
and = everywhere else. In terms of a blast wave model, 
high 71 can come about in highly energetic, fast shocks, 
where a high fraction of the swept up energy is used 
to accelerate particles, and the lower-ener gy part of the 
distribution becomes particle- starved (e.g.. lDermer et ail 
120001 iDermer fc Menonl 120091 ). The hard X-ray spectra 
in some blazars m ay indicate very ha rd electron spectra 
at lower energies (jSikora et al.l 12009). Although blazars 
can be quite variab le on timescales as sh ort as hours (e.g., 
Abdo et al.l I2009H ) or even minutes aronian et al.l 



20071 ). we will assume their average or quiescent emis- 



sion can be described by a steady state solution to the 
electron continuity equation, where continuous injection 
is balanced by cooling and escape. The power continu- 
ously injected in electrons is given by 



(2) 



Linj,e = m e C / C?7 7 Q e {l) 



or, using Equation (J]), 



(7i 2 ^-72~ 9 )/(<Z-2) 
111(72/71) q = 2 



(3) 



We assume an energy-independent escape timescale given 
by 



te 



(4) 



where R' b is the comoving radius of the blob and 77 is a 
constant > 1. In this case, where 71 < 7 C (the slow- 
cooling regime) the electron distribution can be approxi- 
mated as 



iVe(7)~Q0*e S c7c 1 ( 7/7o )-«-l 7C < 7 < 72 • (5) 

If 7 C < 71, i.e., the fast-cooling regime, 

M M~ D i -v / (7/7i)~ 2 7c < 7 < 7i 

N e ( 7 ) ~ QoW7c7i I ( 7 / 7l )-«-i 71 < 7 < 72 ' 

(6) 

Here we assume the electrons are cooled by synchrotron 
emission and Thomson scattering, so that 



3m P c 2 



7c 



4cctt(m' b 



L sy,tot 



+ r 2 



(7) 



is the cooling electron Lorentz factor, where primes 
denote quantities in the comoving frame of the blob. In 
this equation, the magnetic field energy density in the 
blob comoving frame is 



B 2 
8?r 



(8) 



Blazar Sequence 



7 



the total synchrotron energy density is 



, a T u B 

L sy,tot nR ,2 



dj 7 2 We (7) . 



(9) 



(e.g.. IBottcher et aLlll997|) and the external energy den- 
sity u ex t is assumed to be isotropic in the proper frame 
of the AGN. As discussed in Section Q] the exact nature 
of the external radiation field is not known, and may not 
even be the same for all blazars. 

The two main differe n ces b etween our model and that 
of IBottcher fc Dermerl pOOl arc (1) we relax the as- 
sumption that Sd = T; (2) we do not assume that the 
magnetic field energy density is a constant fraction of the 
electron energy density, 



4nR' b 3 /3 



dj 7 N e (j) . 



(10) 



Since spectral modeling of blazars finds that the ratio 
u' B /u' e vary by » 5 orders of mag nitude from source to 
source (e.g.. iGhisellini et aLl feOlO). we do not feel this 
assumption is well-justified. We assume that a fraction 
(r) of the accretion dis k makes up the exte rnal radiation 
field (u ex t). Following IBottcher fc Dermerl (|2002f) . we as- 
sume that t on Ld, the accretion disk luminosity, so that 
u ex t oc rLrf oc L 2 ,. From the B landford-Znajek mech- 
anism (jBlandford fc Znaiekl[l97l . one finds the power 
extracted from the black hole rotation to be L B z oc B 2 t , 
where B m is the magnet ic field in the magnetosp here 
near the black hole (e.g. iCavaliere fc D'Elial [2001 . If 
Ld oc Lbz and B oc B m , where B is the magnetic field 
in the primary jet emitting region, then one finds that 
u ex t oc B 4 . Therefore, we decrease the magnetic field 
assuming 



B = B Q 



( Uext 
V u 



1/4 



(11) 



consistent with IBottcher fc Dermerl (|20Q2| ) . 

In the slow-cooling regime, Equation ([9]) can be inte- 
grated to give 



l sy,tot 



7, 



i f c „ 71 — ^ 

3-q 

qr+1 J2 7c 

2-q 

UTUgQotesc ll~ q 



(12) 



nR'i 



3-q ' 



while in the fast-cooling regime, 

, _ T T u' B Q G teseleli q ~ 



l sy,tot 



nR' 2 



7i 2 (7i"7c)+ (13) 



-q+2 -q+2-\ 
PTu' B QQt e sclc 2-q 

nR'i 71 



For a given set of parameters, the nonlinear nature 
of u' sy means that the above equations do not have 



a simple closed form solution for N e (^). The nonlin- 
ear effects of SSC coolin g have been expl ored in de- 
tail bvlSchlickeiserl (l2fJ09l ); lSchlickeiser et all (|2010D ; and 
iZacharias fe Schlickeiserf (|2010L 12012^ We solve for 
iV e (7) numerically. 

The observed isotropic synchrotron luminosity in 
the frame of the AGN in the (5-approximation ( e.g., 
iDermer fc Schlickeiserll2002t IDermer fc M enon 200J) is 



2(5 



= eL sy (e) = -D C( TTUB7 3 ^e(7) , 



where 



7 



(14) 



(15) 



e is the dimensionless emitted photon energy, cb = 
B/B cr , and B cr = 4.414 x 10 13 G. For 7 < mm(7i,7 c ), 
L" v oc e 4 / 3 , assuming it remains optically thin. The syn- 
chrotron peak will occur at 7 = j c in the slow-cooling 
regime, and at 7 = 71 in the fast cooling regime, so the 
peak synchrotron luminosity will be 



L tl = L Z* = -f^u' B Q,tesc I 



This peak will occur at 



7c 9 7i < 7c 
2-9 

7c7i 7c < 7i 



(16) 



sy 
^pk 



sy 
pk 



m e c 2 ~ 1.23 x 10 20 Hz _ 1)68 \ 7i 7c < 7i 



7c 7i < 7c 



( 17 ) 

The synchrotron luminosity in the frame of the blob, 
L'u is given by Equation (fT6|) with Sx> = 1. From this 



pk 

one can find the peak synchrotron energy density, 

, __ R' b L pk _ u' b <7t n / ll' q 7i<7c 
Usy ' pk ~ c 47ri?; 3 /3 " 2ttR' 2 Wesc 1 7c7i 2 ~ 9 7c < 7i 

(18) 

The SSC luminosity at t he peak in the Tho mson regime 
can be approximated by (jFinke et al.H2008T) 



- s ' s '< _ "^D „„ 



f -..-.I _-L>„,< /, , 1 7c q 7i < 7c 

Pk ~ —CVTUsy^Qotesc <j ^2-q % < ^ , 



or, using Equation ([18 



(19) 



tSSC _ -_2-./ (n t N2 / 7c 2q 7l<7c 



3ttR' 2 



( 2 °) 

The peak luminosity from Thomson-scattering an ex- 
terna l isotropic radiation field in the ^-approximation is 
(e.g.. IDermer fc Schlickeiserll2002h 



TEG_f&„„ „ n + / 7c q 7l<7c / 91 s 

L vk — °D c<J TU e xtQotesc S 2-q „ . (ZlJ 

P I 7c7l 7c < 71 

The Compton dominance (Ac) is given by the ratio of 
the peak Compton-scattered component to the peak of 
the synchrotron component, 



max[L^ c , L^ c ] max[^« el( , u' sy pk ] 



■(22) 
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where we have ignored a bolometric correction term ~ 1. 
Note that SSC emission has the same beaming pattern 
as synchrotron, and thus L^ c /L s p v k does not depend 

on the viewing angle; however, L pk c does not have the 
same beaming pattern as synchrotron and SSC, and so 
L vk/K y k is d e pendent of the viewing angle through <5d 
(|DermeiHl995tlGeorganopoulos et al.ll20Qll ). 

3.2. Results 

We would like to use the model described in Section [3~Tl 
to phenomenologically reproduce the blazar properties 
shown in Figures [2] and [5] Blazars are necessarily not 
observed at the same angle, and we would like to take 
this into account. By contrast, in previous attempts to 
explain the "blazar sequence" (e.g., IGhisellini et alj!998t 
IBottcher fc Dermeril2002T: iGhisellini et al.ll2010ll . authors 
usually assume all objects are viewed at an angle such 
that Sd = r. Further, we would like to do this by varying 
the least number of parameters possible between objects. 
Finally, we note that we attempt to reproduce Figures 
[21 and [5] directly, something previous authors have not 
done. We are attempting to reproduce the population as 
a whole, and not individual blazars. 

Our modeling results are shown in Figures |6] and [Jj 
The curves show the sources at at constant angle, ac- 
cording to the model, so that along the curve only B 
and u ex t vary according to Equation [TT] The model pa- 
rameters are sh own in Table [H and a re fai rly close to the 
values used by IBottcher fe Dermerl ()2002h . The sources 
with low B and low u ext have high z/^, due to the small 
amount of cooling. As B and u ext increase, the cooling 
increases, and hence 7 C decreases. Since u pk oc e p k oc j^B 
(Equation [17]), and 7 C decreases faster than B increases, 
VpL wm decrease. If the cooling is great enough, and 
7c < 7i, the curves enter the fast cooling regime, and 
v' 



the j et of a single source may have velocity gradients 



is associated with 71 instead of 7 C (again Equation 
[TT]). The transition between the fast and slow cooling 
regimes leads to the sharp break at < 10 13 Hz (the ex- 
act location depends on the angle) seen in the curves in 
Figures [5] and [7] In Figure [JJ there is also a break in the 
curves between w 10 14 and ~ 10 15 Hz caused by the tran- 
sition between SSC and EC. Note that for SSC, Ac does 
not depend on angle, while for EC, it does. The curves 
reproduce almost all of the objects in Figures [6] and UJ 
although there are objects with Ac < 0.1 and v s pk > 10 15 
Hz that are not reproduced. A possible explanation is if 
there exist structures within the jet, so that as one views 
a jet more and more off axis, one views a different region 
of the jet (|Mever et al.ll2011f) . This would only work if 
B, Qo, and/or R b also varied not just 8d, since Ac is 
independent of Sd in the case of SSC. 

According to this simple model, all sources have T = 
30. Of course, in reality, sources will not have jets moving 
at the same speed, as s hown by very long baseline inter- 
ferometry (VLBI; e.g. 
20091 IPiner etall 2008 



Jorstad et al.l 120051 : iLister et al.l 
201 Of ), which also typically show 



slower values for T. However, the main emitting regions 
could be on smaller size scales than is possible to re- 
solve with VLBI, and so may ha ve different speed s, al- 
though note that much lower (e.g., Piner et aT1l2008D and 



much higher (e.g.. iMarschc r et al.l 



201(1 " Lorentz factors 



e.g.. iChiaberge et al.l l2000l iStawarz fc Ostrowski 



Georg anopoul os fc Kaza nas 20031: IGhisellini et al.l 12005 
Meyer et al.l 12011^ However, we assume a gradient in 



2002 



r is not necessary to reproduce the general trends (cf. 
IMever et al.l 120 1 If) . and this value may represent an av- 
erage value for the main emitting region. The magnetic 
field strength values, which span from B — 0.050 G to 
B = 8.6 G, are consistent with those found from s pectra l 
modeling of blazars (e.g., IGhisellini et alj 119981 12010(1 . 
These values are also consistent in that the modeling 
generally shows higher B for FRSQs than BL Lacs. 

Blazar SED modeling also indicates that the 7-ray 
emission from FSRQs is likely from Compton scatter- 
ing of an external radiation source, while for HSP BL 
Lacs SSC is a ble to provide a good fit to the 7-ray 
emission (e.g., IGhisellini et"aT][T99l [20Toh . A correla- 
tion between Ac and core dominance, a proxy for 9, 
found in FSRQs, is also evidence that EC dominates 
in these sources ()Mever et al.l [2012'). BL Lacs tend to 
have less prominent "blue bumps" from accretion disks 
and weaker broad lines, two of the leading contenders 
for the seed photon source for external Compton scat- 
tering. The presumed parent population of BL Lacs, 
FR I radio galaxies, als o seem less likely to have dust 
tori (|Donato et al.l [200^1 . another leading contender for 
the external photon source. Wide Field Infrared Survey 
Explorer observati ons of BL Lacs do not show any evi- 
dence for dust tori (|Plotkin et al.|[2012T l . So whatever the 
source is, it seems reasonable to assume it is greater for 
m ore FSRQs, and weaker for BL Lacs. This is similar 
to IBottcher fc Dermerl (|2002l ) but in contrast to Ghis- 
ellini et al., who use a "binary" u ex t where the exter- 
nal radiation field is either "on" or "off" above and be- 
low a certain accretion rate (e.g., Ghisclli ni fc Tavecchiol 
I2008al) . They justify this from reverberation mapping 
campaigns, which show the disk luminosity is propor- 
tional to_the_s^uaje_of the BLR distance from the disk 
(e.g.. iBentz et al. 2006), which would result in the same 
energy density observed, as long as the emitting region 
is inside the BLR, and the same fraction of disk emis- 
sion is reprocessed by the BLR. However, there is no 
guarantee that this fraction is the same in all sources, 
or that if the external radiation source is a dust torus 
that it follows the same relation as the BLR. It seems 
more likely that u ex t can have a range of values, rather 
than just one. Whether the external photon source 
is a dust torus or BLR, its luminosity should be less 
than the Eddington limit for a black hole with mass 
M = 1O 9 M , or L Edd = 1.3 x 10 47 Af 9 erg s" 1 where 
M g = M/(10 9 M Q ). If the radius of the external source 
is Rext = 10 18 cm w 1/3 pc, then the maximum energy 
density will be 



Uext 



->Edd 



0.35 erg cm" 3 M g R 



IS 



(23) 



have been observed. It has also been suggested that 



where = R ex t/(10 18 cm). It is unlikely though that 
the external photon source will be radiating all of the 
black hole's accretion luminosity, so u ext should be lower 
than this by at least a factor 10. We have thus chosen our 
parameters so that the maximum external energy density 
is u ext = 0.035 erg cm~ 3 . 
We have also computed the jet powers in Poynting flux 



Blazar Sequence 



9 



10 

io 4 

l ° 4 

' 5/3 in 4 ' 
M 10 

O 

>» 4 

* io 
io 4 
io 4 
io 4 




• » pi » ool, 

o I 4 



o. 



10 



10 ' 



10 



10 

V 



10" 
'[Hz] 



10 



10 



10 



Figure 6. Same as Figure[2] with curves showing our model plot- 
ted at various angles, shown in the legend. Along the curves only 
B and u ex t are varied. Model parameters are found in Table [3] 
Diamond symbols show the location of the model SEDs plotted in 
Figures \E\ and [9] with the numbers indicating the location of the 
SEDs in Figure E] 




Figure 7. Same as Figure[5]with curves showing our model plot- 
ted at various angles, shown in the legend. Along the curves only 
B and u ex t are varied. Model parameters are found in Table [3] 
Diamond symbols show the location of the model SEDs plotted in 
Figures [8] and [9] 



(24) 



and electrons, respectively, by 

Pj,B = 2TrR' b r 2 c(3u B 
and 

P j>e = 2TrR' b r 2 cPu' e (25) 

(ICelotti fc Fabian|[l99l ICelotti et all [20071: iFinke et al.l 
I2008H . These results are also shown in Table 03 The 
jet powers are consistent with those fou nd previously by 
other autho rs for spectral modeling (e.g. JGhisehini et"aH 
[2OO9ll2OT0h . For sources with high B, naturally Pj t B is 
higher. In the slow cooling regime, Pj te is nearly inde- 
pendent of B and u ex t , decreasing with increasing B and 
u ex t only slightly. In the fast cooling regime (Equation 
[6]), N e (pi), and hence u e , is directly proportional to j c , 
so that as B and u ext increase, j c decreases, leading to 
lower u e . 

The SEDs of several blazars along the sequence, seen at 
an angle to the line of sight 9 — and 9 — 2° can be seen 



Table 2 

Parameters (described in the text) for 
theoretical curves in Figures [6] and [7] 



Parameter 


Value 


71 


10 2 


72 


10 6 


Rb 


10 16 cm 


V 


1.0 




2.0 X 10 42 erg s" 1 


Qo 


6 X 10 48 s" 1 


1 


2.5 


u 


4.0 X 10~ n erg cm -3 


r 


30 


Bo 


0.05 G 



in Figures [5] and [51 respectively. These SEDs correspond 
to the diamonds in Figures [5] and [7J The parameters for 
these curves can be found in Table [3l The external ra- 
diation field is assumed to be isotropic in the galaxy's 
frame and monochromatic, with dimensionless seed pho- 
ton energy e = 5 x 10~ 7 in the galaxy's frame, which 
is around what one would expect for the peak emission 
from a dust torus with temperature Td us t ~ 10 3 K. The 
calculations were performed using the exact synchrotron 
emissivity and the full Compton cross section, accurate 
in the Thomson through Klein-Nishina regi mes. The de- 
tails o f th ese calculations can be found in IFinke et al.l 
( 20081) and lDermer et~ail ((1)09). Since exact expressions 
are used, the parameters L S X, i/fj£, and Ac differ slightly 
from the values found in Figures [S] and [7] Note that 
synchrotron-self absorption is included in these curves as 
well, and that it is more apparent for the higher power 
blazars, where the magnetic field is larger and the emit- 
ting region becomes self-absorbed at higher frequencies. 
A contribution from an underlying accretion disk is not 
included, although this has been shown to dominate the 
optical continuum for FSRQs and some BL Lacs. These 
curves ap pear similar to SEDs observed from blazars. 
Similar to Bottchcr et al. (2002), we do not attempt to 
reproduce any individual blazar. However, to demon- 
strate the similarity to actual blazar SEDs, we plot data 
for several blazars on these Fig ures as well. The s e inclu de 
the CRATES 0630-2406 from iPadovani et all (120121) in 
Figure and the Mrk 421 S EP fromlAbdo et all (1201111 
and the low state 3C 279 from lHayashida et al.l ( 20121 ). in 
Figure [9l The symbols have the same color as the curve 
which is the closest match. The curves are not a per- 
fect fit to the SED data, although they are a reasonable 
representation. 

The scenario described here predicts that a large frac- 
tion of FSRQs are emitting in the fast-cooling regime. 
In this regime (j c < 71), one expects the electron in- 
dex below the peak (associated with 71) to be p = 2. 
as shown in Equation ©. For FSRQs, Compton scat- 
tering by these electrons should make hard X-rays (~ 
10 - 100 keV) and the spectral index from Thomson 
scattering would then be Ti — (p + l)/2 = 1.5, assum- 
ing one is observing reasonably far below the break. In 
the Swift-HAT spe ctra from 36 months of observations 
(|Aiello et al.ll2009f ), almost every FSRQ has an X-ray 
spectral index consistent with Tj = 1.5 within the er- 
ror bars, and the mean for t he sample is 1.6 ± 0. 3. The 
BeppoSAX six year catalog (|Donato et al. I [20051 ) has a 
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Figure 8. SEDs for the sequence seen in Figures [6] and [7] as ob- 
served at 6 = 0. The numbers correspond to the numbered di- 
amonds in Figu re [6] The circle s are t he SED data for CRATES 
0630-2406 from Padovani et al.l l laHlT ). 




v [Hz] 

Figure 9. SEDs for the sequence (diamonds) seen in Figures [6] 
and [7] as obser ved at 6 = 2°. T he circles are the SED data for 
Mrk 421 from lAbdo et al.l H201U). and the squares are the SED 
data for 3C 279 from Havashid a et al.l (USD- 
large number of sources are near I\ w 1.5 with a mean 
spectral index of 1.59 ± 0.05, although not all sources in 
that catalog are consistent with I\ = 1.5. The diversity 
of hard X-ray spectral indices could be explained by the 
location of the Compton component within the BAT or 
BeppoSAX bandpass. Softer spectra (1^ > 1.5) could be 
those with their Compton peak being in or near the BAT 
or BeppoSAX bandpass, or by a contribution from SSC 
emission. Harder spectra (r, < 1.5) could be caused by 
the photons scattered by electrons below j c being in the 
hard X-ray bandpass. 

4. SUMMARY AND DISCUSSION 

In Section [5] we have found anti-correlations in the 
2LAC that were found previously by many authors in 
other samples. The plot and correlation with Ac in par- 
ticular is novel, since the 2LAC allows for a larger sam- 
ple of the Compton component than any previous cat- 
alog. We include on this plot for the first time a large 
number of blazars without known z, and showed that 
the anti-correlation with v pk is not a relic of ignoring 
blazars which lack spectroscopic redshifts. This relation- 
ship see ms to have a physic al ori gin. 

Both iMever et all (|2012D and IGhisellini et al.l (|2012H 



plot the LAT 7-ray luminosity (L 7 ) versus the LAT spec- 
tral index (F 7 ; a proxy for v^ k ), and showed that the four 
sources from iPadovani et "all (|2012f) with high v s p \, high 
L s pk do not have excessively large L 7 or excessively small 
T 7 , and are consistent with other LAT blazars. However, 
it is possible that the redshifts of many BL Lacs in the 
future could be measured, and they could be found with 
high L 7 and low F 7 . So it is not possible to assess this 
selection effect with a plot of L 7 versus L 7 . 

In Section|3]we developed a simple model to explain the 
plots and anti-correlations from Section [21 This shape is 
explained by the increased external radiation field energy 



sy 



is asso- 



density and magnetic field strength, where v pk 
ciated with either "f c (the Lorentz factor of the break in 
the electron spectrum due to cooling) for sources in the 
slow-cooling regime, and 71 (the lowest Lorentz factor 
of the injected electrons) for sources in the fast-cooling 
regime. 

This theory is quite simple, and neglects numerous 
important effects. Perhaps most importantly, blazars 
are highly variable, and we treat them with a steady- 
state solution to the electron continuity equation. In 
any study of a large population of blazars, however, this 
is almost unavoidable. In developing this theory, ap- 
proximate expressions were used, neglecting the exact 
synchrotron emissivity and Compton cross section, es- 
pecially Klein-Nishina effects. However, a comparison 
between the predictions and exact calculations (Figure 
[5]) shows that these approximations seem reasonable. 

Despite its simplicity, this model makes a number of 
predictions. It predicts that the single most important 
parameter in determining the luminosity (L s p y k ) of blazars 
is 9, mostly independent of the frequency of their syn- 
chrotron peak (Figure [S]). This also means that as fainter 
sources are found, these sources will not have excessively 
large values for Ac, since large values for this quantity 
are due to EC, and Ac is strongly dependent on 9 if EC 
is dominant. There should not be any sources found with 
L s p v k < 10 45 erg s^ 1 and Ac > a few. This model also 
predicts that no "blue FSRQs", that is, no HSPs with 
significant BLR or dust torus luminosity will be found. 
It predicts that any source with high Ac (> a few) will 
be an LSP and it will be in the fast cooling regime, so 
that the spectral index of the EC emission below the 7- 
ray peak will be Ti ~ 1.5. For most sources, this will be 
in the soft 7-rays (< 10 MeV) down to the hard X-ray 
regime, perhaps as low as 10 keV in some cases. 

This model predicts that if the BL Lacs sources be- 
come more redshift-complete, more sources with high L s pk 

and v s pk ("bright HSPs") will be found, although these 
sources will not be as bright as the brightest LSPs (Fig- 
ure IS]) . These sources will have their jets highly aligned 
with our line of sight. Note that our scenario differs 
somewhat from the one of IGhisellini eFall (|20Tl . They 
also predict that bright HSP sources will be found, but 
in their case, these sources are "blue FSRQs", with the 
primary emitting region found outside the BLR, avoid- 
ing an EC component and a large Ac, since the 7-ray 
emission will be due only to SSC. Their scenario and 
ours would produce essentially identical SEDs for bright 
HSPs, so distinguishing them is not possible based on 
SEDs alone. If it could be found that these sources have 
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Table 3 

Parameters (described in the text) for curves shown in Figures [8] and [9] 



Parameter 1 2 3 4 5 6 7 

B [G] 0.050 0.12 0.28 0.67 1.6 3.8 8.6 

u ext [erg cm" 3 ] 4 x lO" 11 1.3 x 10" 9 4.1 x 10~ 8 1.3 x 10~ 6 4.2 x 10" 5 1.4 x 10~ 3 3.5 x 10~ 2 

7 C 1.4 x 10 5 3.9 x 10 4 1.1 x 10 4 2.7 x 10 3 5.4 x 10 2 50 2.7 

Pj,B [ergs" 1 ] 1.7 xlO 42 9.6 x 10 42 5.4 x 10 43 3.1 x 10 44 1.7 x 10 45 9.8 x 10 45 5.0 x 10 46 

P 3 ' e [erg s" 1 ] 1.3 x 10 45 1.3 x 10 45 1.2 x 10 45 1.2 x 10 45 9.5 x 10 44 4.5 x 10 44 7.6 x 10 43 



their jets highly aligned, that they do not have signif- 
icant BLRs or dust tori, or that the BLRs are not the 
seed photon sources for EC in "red" (i.e., LSP) FSRQs, it 
would favor our scenario. All of these things are difficult 
to determine observationally, however. The nonthermal 
synchrotron makes broad emission lines or a dust compo- 
nent nearly impossible to observe, although jet alignment 
may be possible to dete rmine with VLBI observations 
(e.g., Uorstad et aHl2005D . 

Our model c an also be contr a sted wi th the Monte Carlo 
simulations of iGiommi et al.l j2012a). In their simula- 
tions, they randomly draw for each blazar several prop- 
erties including 7 p fc (the peak of the electron distribution) 
Sjj, and the strength of the broad emission lines, all in- 
dependent of each other. Although they do not explore 
the Compton dominance in their paper, if these proper- 
ties are in reality independent of each other, and if the 
BLR is the se ed photon sour c e for C ompton scattering, 
the scenario of lGiommi et al.l ((2012a) would not produce 
a correlation between Ac and u^, and would produce 

objects with high Ac and high z/^. Since no such ob- 
jects are found (Figured]) and a correlation between Ac 
and v s ^ k is found (Table [T]), their scenario does not seem 
to be in agreement with observations. 

Our simple model assumes the primary emitting 
region has a single Lorentz factor. Ba sed on op- 
tical ( Chiaber ge et al.l 120001 ) and 7-ray (jAbdo et al.l 
l2010bl lal) observations of radio galaxies, it seems highly 
likely that jets are stratified in speed either per- 
pend icular (jStawarz fe Ostrowskil 120021 : iGhisellini et al.1 
12005ft or parallel to the jet's direction of motion 
( Georganopoulos fe Kazanasll2003l ). Indeed. IMever et all 
(2011) have an explanation for an "L" shape in a plot of 
Lpi versus v"%, based on sources with varying Lorentz 
factors within a source. In our study we neglect radio 
galaxies and stratified jets. However, modeling of ra- 
dio galaxies (iChiaberge et al.1 120011 : lAbdo et all l2009cL 
I2010ct iMigliori et al.ll2011h indicates that these types of 
structures almost certainly exist. Whether they can ex- 
plain the properties of blazars as well as radio galaxies 
remains to be seen. 

Another possible way to decipher the correct model 
could involve estimating the power injected into the jet. 
This power can be related to the power needed to cre- 
ate a cavity in the hot X-ray emitting I CM surrounding 
radio galaxies ([Birzan et al.ll2004l 120081 : iCavagnolo et al.l 
l2010l) . This power seems to b e correlated with th e ex- 
tended lobe's radio power, and IMever et al.l ()2011f ) have 
used the radio power as a proxy for jet kinetic power. 
This angle-independent measure of the jet kinetic power 
may provide a way of distinguishing viewing angle effects 
from intrinsic power effects. 



As a redshift-independent quantity, the Compton dom- 
inance is a useful tool for exploring blazar properties, 
including the large number of BL Lacs without known 
redshifts. The large new blazar catalog from the Fermi- 
LAT is able to characterize the Compton component for 
a larger number of objects than previously possible, mak- 
ing it valuable for determining its relationship to the 
blazar sequence. 

We are grateful to CD. Dermer, M. Georganopou- 
los, M. Lister, and K. Wood for useful discussions on 
the blazar sequence, to CD. Dermer and A. Wehrle 
who gave helpful comments on a previous version of this 
manuscript, to M. Hayashida for the SED of 3C 279, 
and to the anonymous referee for very useful comments. 
This work was partially supported by Fermi CI Grant 
NNH09ZDA001N. 
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